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^ ' Abstract 

o : 1 - 

Q^ . We present a calculation of the mass of the Sq pseudoscalar be [Be) state us- 

^» I ing a non-perturbative measurement from quenched lattice QCD. We find Mbc = 

"^ ■ 6.386(9) (98) (15) GeV where the first error is statistical, the second systematic due 

'~j^ . to the quark mass ambiguities and quenching and the third the systematic error due 

^Li, to the estimation of mass of the %. 
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X ■ 1 Introduction 

Little is known experimentally of the properties of the be spectrum of states, apart from 
the detection of the pseudoscalar ^5*0, be (Be) state by the CDF collaboration |1[]. At the 
very least these states provide an excellent "blind" test for some of the techniques used 
in measuring non-perturbative properties such as the mass spectrum and decay widths. 
The study of non-degenerate heavy quarks provides us with a unique probe of QCD. The 
weak decays of be states provide a plethora of new methods for calculating the less well 
determined CKM matrix elements, in particular its semi-leptonic decay mode 

B.^J/tfjlV , (1) 

will provide an accurate measurement of |\4fe| [0; 0|- 

The spectrum of states with the ground state mass subtracted, when calculated via 
potential models 0, |^, for example, is in broad agreement with quenched lattice QCD 
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0, 0, . Despite concerns that potential models would have large relativistic corrections (as 
a simple centre-of-mass argument implies (f ^/c^) for the charm quark is approximately 0.5) 
it appears that it is still a reasonable method for calculating the spectrum. Nonetheless, a 
non-perturbative measurement of the spectrum is important and is being actively pursued. 
The potential model approach unfortunately does not provide, a priori, a method for 
determining the ground state mass and hence a number of different phenomenological 
methods have to be used. Lattice QCD, on the other hand, can provide an ab initio 
method for determining its mass. 

The use of a non-relativistic hamiltonian to describe the behaviour of a bottom quark 
in quenched lattice QCD has been employed very effectively for the spin-independent 
sector of the bb spectrum [|l^ . For relatively coarse lattice spacings, this technique is also 
applicable for cc states and indeed two calculations [^ ^^ of the be spectrum and rriBc 
using non-relativistic b and c quarks have been carried out. However, as the cc spectrum 
is not reproduced very accurately using this approach []T2[ there remains some question as 



to its accuracy for the be sector. Another approach which will work for lattice spacings 
where aruc < 1 is to use a discretised relativistic Lagrangian, which, for the lattice spacing 



used in this letter, reproduces the charmonium spectrum reasonably well [|1^, |TJ]. Here we 
take the approach of combining these separate methods for treating the b and c quarks. In 
essence this involves computing Green's functions from a discretised NRQCD hamiltonian ( 
NRGF ) with propagators calculated using a improved discretised QCD Lagrangian devised 
originally by Sheikholeslami and Wohlert |T^ ( SWP ) from the same set of quenched ( zero 
flavour ) QCD gauge configurations and then combining them into two point functions. 
The damped exponential behaviour of this correlation function will provide information 
about the mass of the states we are interested in. 

The rest of this letter is organised as follows : in the following section we describe 
a method for computing msc directly from zero momentum correlation functions, and 
demonstrate why this has the smallest systematic error; after this the computational details 
of the calculation are outlined, explaining briefly the details of NRQCD Hamiltonian used 
for the bottom quark and the Lagrangian used for the charm quark; we then present the 
numerical results and in the conclusions we compare the results with other calculations of 
mB.. 

2 Ground state mass definitions 

For large Euclidean times, the two point function of a meson operator constructed from 
either NRGF's or SWP's will take the following form 

d^x(0|M(x,t)Mt(0)|0) oM/mi Aexp(-amit/a) , (2) 

where the integral sign can also indicate a sum over discrete lattice points. (We have 
also discarded the possible effect of the finite size of the lattice in the time direction 



for relativistic correlation functions.) As demonstrated in |]T8[ and |]T9[ the variable ami 



(where a is the lattice spacing), also referred to as "the pole mass", has quite different 
interpretations depending on whether the operator M(x, t) is constructed from NRGF's 
or SWP's. The former is the binding energy of the state with the heavy quark mass 
subtracted; the latter satisfies the relationship 

arrii^^ = aM + ac{amQ,a,V) + ad{a,V) , (3) 

where M is the mass of the state, ad are the systematic errors from the finite lattice spacing 
and volume and ac is a lattice artefact due to the breaking of Lorentz invariance on the 
lattice This parameter vanishes for amq << 1 and grows to be a large correction as airtq 
approaches 1 and beyond. Nonetheless, the difference of the pole masses, for either NRGF's 
or SWP's will reproduce the physical mass difference (assuming other discretisation effects 
are under control). Since this process of changing the zero of the quark mass is additive, 
one therefore expects 

mis. - l/2(mf^f^^ + ml^/) = Mb. - 1/2(M^, + M,J , (4) 

where the RHS is the physical mass difference and miBc is computed using correlation 
functions with the composition used in this letter. As one can easily construct a B along 
the same lines one can also state that 

miBc - (miB + mf^^) = Mb. - {Mb + Mj,) . (5) 

It is therefore possible derive an expression for Mb. simply from the pole masses without 
resorting to using other possible mass definitions. 

The mass of the rji, has yet to be determined, and has to be estimated from the T and 
a potential model calculation of the hyperfine splitting. A larger error is that resulting 
from the quenched approximation. A first estimate of the effect of quenching is from 
the variation of the effective lattice spacing between observables having different "typical" 
momentum scales. In this light we estimate this from the difference between the masses 
determined from equations § and ^ 

A more traditional approach to calculate the mass of the Be is to determine the "kinetic" 
mass, (often labelled am2) which is derived from the dispersion relationship of pole masses 
at different finite momenta. This however has the drawback that it has a much larger 
statistical error than the pole mass. In particular, while the kinetic mass is proportional 
to the inverse of a difference (i.e. the zero and smallest non-zero momentum pole masses), 
the mass derived from equations @) and (^) are directly proportional to it. 

3 Computational Details 

The calculation was carried out in quenched lattice QCD using the Wilson gluon action at 
(3 = 6.2 on a 24^ x 48 lattice. Technical details on how these configurations were calculated 



and the exact form of the gauge action can be found in [^ . In terms of the string tension 



resulting from the static quark potential, the lattice spacing is approximately 0.07 fm, 
resulting in a spatial side length of 1.68 fm. 

The non-relativistic Hamiltonian used was corrected to 0{mv'^). The coefficients in the 
Hamiltonian were tree-level mean-field estimates. The mean field improvement coefficient 
is calculated from the plaquette. Further details on the Hamiltonian used are described in 
0]. The bare mass parameter aM^ for the he states was chosen to be 1.22. From ||T0| this 



choice of the bare quark mass parameter is consistent, within quenching errors, with the 
b-quark mass when defined by the kinetic mass of the T, and the scale determined from 
the T' — T splitting. As the splitting of equation (H) should be most sensitive to a variation 
in the 6-quark mass, we use this scale to also fix the charm quark mass. In the case of the 
5-like states we chose a range of values for aM^ from 1.1 to 1.3. A linear interpolation is 
carried out to determine the pole mass for the above bare parameter. 



The relativistic Lagrangian used was the Sheikholeslami-Wohlert action |Tj]. For the 
D, cc and be states, the value of the improvement coefficient, csw, in this action was de- 
termined by the tadpole improvement procedure [^. For the B states, a slightly different 



choice of csw was used, determined by imposing PCAC for light quark masses [^. For 
this very fine lattice spacing, the masses derived from these different prescriptions for csw 



are indistinguishable |23]. For the B and D, a range of light quark masses were used, and 



then extrapolated to the chiral limit. 

In order to determine the charm quark mass we use the kinetic mass of the Dg ^^, [14 
(requoted in table ^ for clarity). This is obtained from correlation functions with one SWP 
with either k = 0.126 or 0.132 and one SWP with a n corresponding to the strange quark 
mass. It can be demonstrated that the kinetic masses of these correlation functions differ 
from the pole masses by only a few percent, so ambiguities in the definition of the charm 
quark mass are minimised. In particular we find that the charm quark mass is close to a 
bare quark mass equivalent to k = 0.126, so the interpolation is relatively mild. 

In the case of the 6c-like states. Coulomb gauge fixed correlation functions were com- 
puted using sources optimised to select the ground ^5*0 state and its first radial excitation. 
The typical radius for these functions lay in the range of 3 to 5. For k = 0.126 the correla- 
tion functions were computed with 201 configurations and for k = 0.132 140 configurations 
were used. As a Green's function calculated using the NRQCD Hamiltonian only contains 
upper spin components, at k = 0.126 the gauge field was also time reversed so as to obtain 
correlation functions using the lower spin components of the relativistic propagator and 
improve the statistics. 

Correlation functions for other states in the spectrum, that is the '^5'i, ^Pi and '^P(2,i,o) 
were also calculated. The absolute binding energies of the ^Pi and ^-Pi,o were calculated 
while the ^5*1 and ^P2 energies were determined from single exponential ratio fits relative 
to the ^5*0 and ^Pq respectively. 

For the B-\ike states, gauge invariant correlation functions were computed with 68 
configurations. Smeared sources were computed using Jacobi smearing |]16| . 



The fit criteria for obtaining the pole masses of be states were the same as those as used 
in p2[, using multi-exponential fits. The 5-like states were fitted to a single exponential. 



The gauge configurations and SWP's were calculated on a Cray T3D at the Edinburgh Par- 



allel Computing Centre (EPCC) while the Coulomb gauge fixing, NRGF's and correlation 
functions were computed on a Cray J90 at the same site. 

4 Results 

4.1 The ground state mass 

The pole masses calculated for the various heavy-heavy states are listed in in table |l|, as 
well the pole masses for charmonia listed in |]TB|, |I1|. The "pole mass" for the T obtained 
from degenerate NRGF's is listed in table |l| [|l^. The chirally extrapolated results for D 
and B are quoted in tables ^ and ||. 

From [|l^ the choice of the bare quark mass parameter of 1.22 is consistent, within 
quenching errors, with the b-quark mass. On the other hand, the choice of k = 0.126 for 
the relativistic quark is not quite consistent with the charm quark mass. By performing a 
slight linear extrapolation on the mass difference defined in (^) using the kinetic mass, we 
find 

amiB, - -{amir,, + ami^J = 0.0557(13) , (6) 

and 

aniiB, - {aniiB + aniiD) = -0.2480(74) . (7) 

As we can see, these two differences vary by a large degree, having a different sign from 
each other. When converted into a dimensionful scale, they vary by approximately 1 GeV. 
In the case of equation (^), the choice of scale is reasonably clear; in order to be consistent 
in the determination of the quark masses, one should employ the scale defined via the 
T' — T splitting, which for this choice of action and parameters is a~^ = 3.52(14) GeV JIU|. 



In the case of equation (|^), the choice is less clear. For the heavy-light masses one would 
normally use a different definition of the lattice spacing, for example that defined from 
the string tension (ideally one would use the 25 — 15 or IP — 15* split for the appropriate 
heavy- light combination, but the statistical error is too large in this case). Nonetheless the 
above argument for the heavy-heavy difference indicates that the lattice spacing from the 
T' — T splitting should be used here as well. The resulting estimate for Mb, from equation 
(0) varies by approximately 200 MeV solely by the choice of this ambiguity in the lattice 
spacing. The result for the mass from equation (^ lies at the centre of this interval. We 
therefore take the mass determined from equation (^ as our best estimate of the central 
value and take the variation with respect to the scale for equation (|^) as a conservative 
estimate of the systematic error due to quenching and quark mass ambiguities. 

We also require an estimate for the hyperfine splitting for T — 77^ as experimental data 
for Mrjb is not yet available. Quenched lattice studies of the charmonium spectrum have 
always underestimated the hyperfine splitting (as the central charge from the Coulombic 
term is smaller than usual for QCD where nf = 0). In order to estimate the size of this 



splitting we use the value of the hyperfine splitting determined in [|TU| in conjunction with 



the the potential model results quoted in |^ . From this we estimate the T — rjb splitting to 
be 60(30)MeV. Hence we take M^^ = 9.400(30) GeV. Using these parameters, one finds : 

Ms, = 6.386(9)(98)(15) GeV , (8) 

where the first error is statistical, the second is the systematic error due to quenching 
and the third is the systematic error in the estimation of M^^ 

4.2 The spectrum 

The absolute or relative binding energies for the higher states of the spectrum are listed 
in table |. In figure we plot these results along with the final lattice results at /5 = 5.7 
quoted in P] and the potential model predictions in |^, with the scale determined from 
the spin averaged IP — IS splitting, which has the smallest statistical errors. In general 
these results are in broad agreement with each other, although there does appear to be a 
significant difference in the hyperfine splitting between the lattice results and the potential 
model calculation. 

5 Conclusions 

In this letter we have presented a calculation of the ground state mass of the B^ and the 
immediate spectrum of states. Although the calculation has been carried in the quenched 
approximation, by using two different possible derivations of the ground state mass which 
use a wide range of different possible quark masses, we believe we have conservatively 
estimated its effect. 

The results for the spectrum are, for the most part, in broad agreement with the 
results of |^ and [Q. The ground state mass quoted here is consistent with previous 



lattice calculations ||g, |TT|, |2^ in the quenched approximation calculated using the kinetic 
mass definition, using NRGF's for both the charm and bottom quarks on coarser lattices. 
QCD sum rule calculations which are listed in [^ predict that Mbc lies in the interval 
of 6.3-6.5 GeV. Kwong and Rosner |2^ surveyed the phenomenological techniques used 



in determining Mg^ and estimated it to lie in the range of 6.194 GeV to 6.292 GeV. 
Most phenomenological estimates for the mass have been taken to lie in the centre of 
this interval [Q although more recently other work by Fulcher has suggested it to be 

+15 

slightly larger at 6.286 - 6 GeV. Nonetheless, the mass we have determined is certainly 
consistent with the present estimate of the mass from the CDF collaboration |Q of Mbc = 
6.40 ± 0.39{stat.) ± 0.13(sys.)GeV. 

Despite the systematic error due to the hyperfine splitting in bb sector, the largest error 
is the effect of quenching and quark mass ambiguities. Eliminating this approximation by 
the use of dynamical configurations, where sea quark effects are included appears to be 
the most important step in reducing the error on this mass. Likewise, a more complete 
understanding of how the differences of equations (^) and (^) behave as a function of the 
heavy quark masses also needs to be carried out. Such configurations (albeit on coarser 



lattices than used here) are now available and a measurement of the spin-averaged Be 
would be the next course of action. 
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Figure 1: A comparison of the spectrum of states in this letter with other calculations. 
The results denoted by the square and burst were originally quoted in, respectively in 
and 0]. The diamond points represent the results determined in this paper. To eliminate 
the scale for the lattice calculations, the spin averaged IP — IS* splitting has been used. 
The 2E and 2T labels for the ^P2 states represent different possible representations in the 
correlation functions that have overlap with that state. 



K 


aMb 


ami{cc) 


ami{bc) 


ami{bb) 


0.126 


1.22 


1.035(1) 


0.7207(9) 


0.3027(3) 


0.132 


1.22 


0.661(1) 


0.5532(18) 





Table 1: Pole masses determined from correlation functions composed of SWP's or NRGF's, 

for heavy-heavy states. 



K 


ami{D) 


0.132 


0.4332(15) 


0.126 


0.6424(19) 



Table 2: Chirally extrapolated pole masses determined for D-like states. 
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aMb 


ami{B) 


1.1 


0.3192(92) 


1.2 


0.3354(96) 


1.3 


0.347(10) 



Table 3: Chirally extrapolated pole masses determined for 5-like states. 



K 


a(T' — T)/am2(cs) 


0.126 


0.255(16) 


0.132 


0.389(23) 



Table 4: Inverse kinetic masses of heavy-strange correlation functions composed SWF's 
and the heavy is in the region of charm. The scale determined from the T' — T splitting 
is chosen in order to maintain consistency with the b-quark mass definition. The physical 
ratio is 0.283(3). 





K = 0.126 


/t = 0.132 


ami^{m2{Ds)) 


aE(^Si) - aE{'So) 


0.0127±j| 


0.0145±t 


0.0131±^ 


aEii2ySo) 


0.886±^ 


0.756±ii 


0.854±6 


aEi{2fSi) 


0.906±^ 


0.775±9 


0.876±t 


aE{^P^) 


0.864±| 


0.701±^ 


0.822±^ 


aEi'Po) 


0.852±^ 


0.691±^ 


0.814±^ 


aE^P^T) - aEi^Po) 


0.021±2 


0.018±^ 


0.020±3 


aEi^P^E) - aEi'Po) 


0.021±2 


0.027±^ 


0.022±2 


aE{^P,) 


0.861±| 


0.703±^ 


0.822±f 



Table 5: Spectrum of states for bare quark masses and the interpolated estimate at the 
charm mass. 
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